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Abstract. We provide a quantitative picture of non-conserved interface growth from 
a diffusive field making special emphasis on two main issues, the range of validity of 
the effective small-slopes (interfacial) theories and the interplay between the emergence 
of morphologically instabilities in the aggregate dynamics, and its kinetic roughening 
properties. Taking for definiteness electrochemical deposition as our experimental field 
of reference, our theoretical approach makes use of two complementary approaches: 
interfacial effective equations and a phase-field formulation of the electrodeposition 
process. Both descriptions allow us to establish a close quantitative connection between 
theory and experiments. Moreover, we are able to correlate the anomalous scaling 
properties seen in some experiments with the failure of the small slope approximation, 
and to assess the effective re-emergence of standard kinetic roughening properties at 
very long times under appropriate experimental conditions. 
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1. Introduction 

Many natural and technological growth systems exist whose dynamics results from the 
interplay between diffusive transport of aggregating units within dilute phase, and 
their attachment at a finite rate to a condensed growing cluster. Examples abound 
particularly within Materials Science, e.g. in the production of both amorphous (through 
Chemical Vapor Deposition [lj or Electrochemical Deposition [2]), and of epitaxial (by 
Molecular Beam Epitaxy [3]) thin films. Still, a similar variety of surface morphologies 
can be found in diverse systems, such as e.g. bacterial colonies jl], to the extent that 
universal mechanisms are expected to govern the evolution of these non-equilibrium 
systems [5j. Moreover, many times it is not matter, but e.g. heat that is transported, 
but obeying precisely the same constitutive laws (one-sided solidification). 

The objects of interest in the time evolution of these systems are the morphological 
properties of the aggregate which is produced. Thus, a variety of behaviors ensue, 
such as pattern formation and evolution [5J, and surface kinetic roughening [6], whose 
understanding has driven a large effort within the Statistical Mechanics and Nonlinear 
Science communities. In this context, continuum models have played a prominent 
role, due to the compact description they typically provide of otherwise very complex 
phenomena, and their ability to explore generic or universal system properties [7|. For 
the diffusive growth systems mentioned above, the simplest continuum formulation takes 
the form of a moving boundary problem in which particle diffusion in the dilute phase 
is coupled with a boundary condition in which the growth velocity of the aggregate is 
prescribed as, roughly, the fraction of random walkers that finally stick irreversibly to it. 
Thus, the evolution of a growing surface from a vapor or the electrochemical growth of a 
metallic deposit under galvanostatic conditions can be described by the same nonlinear 
system of stochastic differential equations (see Sec. [2] for details). 

Electrochemical deposition (ECD) has been widely studied mainly due to the 
apparent simplicity of the experimental setup and the capability to control different 
relevant physical parameters. Most of the studies in ECD address the growth of quasi- 
two dimensional fractal aggregates. This can be achieved by confining a metallic salt 
between two parallel plates with a small separation between them. Thus, although the 
formulation of the problem can be accurately posed in order to include all the physical 
mechanisms involved, its solution is far from trivial (actually it cannot be obtained even 
in the deterministic case). Moreover, the numerical integration of the corresponding 
equations often involves the use of complex algorithms that take into account the motion 
of the boundary and which, in practice, are difficult to implement. Besides this intrinsic 
complexity, the role of fluctuations and their inclusion in the numerical scheme is still 
an open question. In the present context, the interplay among different transport or 
electrochemical mechanisms is not clear and the predictive power of different models is 
limited to the description of the concentration of cations in the solution or the qualitative 
morphologies of the fractal aggregates [HE]- 

In this scenario, two possibilities arise from a theoretical point of view. On the 
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one hand, one can try to obtain as much information as available from customary 
perturbation methods. Thus, linearization, Green's function projection techniques or 
weakly nonlinear expansions provide important information at both short (linear regime) 
and asymptotic times [10J. Briefly, this approach allows to explore the physics of the 
problem in the so-called small slopes approximation. On the other hand, one can try to 
integrate numerically a more computationally efficient version of the moving boundary 
problem in order to accurately determine the validity of the approximations in the model 
formulation, and also to gain insight into the intermediate time regimes. 

This second approach can be achieved by the introduction of an equivalent phase- 
field model pj] version of the moving boundary problem. Phase-field models have 
been successfully introduced in the last few years to understand moving boundary 
problems which are harder to integrate numerically (or in combination with discrete 
lattice algorithms [12] )• Moreover, they have become a paradigmatic representation of 
multi-phase dynamics whose application ranges from solidification [T3l [Hj or Materials 
Science pH HE1 CE3 EE] to Fluid Mechanics pHED] or even Biology [21]. 

The aim of this paper is to provide a quantitative picture of non-conserved growth 
in which material transport is of a diffusive nature. Working for the case of ECD, we 
will focus on two important and complementary tasks: assessing the validity of the 
effective small-slopes (interfacial) theories, and the emergence of unstable aggregates 
that cannot be described by single valued functions. More importantly, we will show that 
our model is capable to provide both a qualitative and a quantitative picture of growth 
by electrochemical deposition. In particular, our present study allows us to provide more 
specific reasons for the experimental difficulty in observing kinetic roughening properties 
of the Kardar-Parisi-Zhang universality class, in spite of theoretical expectations. 
Finally, we will see that the frequently observed anomalous scaling properties arise 
in this class of systems associated with the failure of the single-valued description of the 
interface, and may yield asymptotically to standard scaling behavior, provided a well 
defined aggregate front reappears. 

Our paper is organized according to the following scheme: in Sec. II we introduce 
our model and discuss in depth the relationship between its defining parameters and 
physical ones. In Sec. Ill, we discuss the validity of the small slope approximation 
and conclude on the need of a full phase-field description of the problem. Section IV 
collects the results of a detailed comparison with specific ECD experiments that allows 
us to derive a consistent picture of pattern formation within our non-conserved growth 
system. Finally, Sec. IV concludes the paper by summarizing our main findings and 
collecting some prospects and open questions. 

2. Continuum model of galvanostatic electrodeposition 

Recently we have proposed a stochastic moving boundary model to describe the 
morphological evolution of a variety of so-called diffusive growth systems [10J. This 
is a large class of systems in which dynamics arises as a competition between supply 
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of material from a dilute phase through diffusion, and aggregation at a growing cluster 
surface through generically finite attachment kinetics. The deterministic limit of the 
model has a long history in the field of thin film growth by Chemical Vapor Deposition 
[221 [231 [21] , where diffusing species are electrically neutral, conceptually simplifying the 
theoretical description. Naturally, in the case of aggregate growth by electrochemical 
deposition (ECD) in galvanostatic (constant current) conditions, the situation is more 
complex due to the existence of two diffusing species (anions and cations) and an 
imposed electric field. Nonetheless, an explicit mapping of the corresponding dynamical 
equations can be made to the CVD system, by means of the following assumptions: (i) 
zero fluid velocity in a thin electrochemical cell; (ii) anion annihilation at the anode, 
whose position is located at infinity, namely, the aggregate height is much smaller than 
the distance separating the electrodes; (iii) electroneutrality away from the electrodes; 
(iv) zero anion flux at the cathode. Within these assumptions, the model ECD system 
takes the form [10] 

d t c = DV 2 c - V • q, (1) 

Dd n c = k D {c-c° eq + TK + x) +q-n, (2) 



V n =Q 



Dd n c - V s ■ J s - q ■ n - V s • p , (3) 



lim c(x, z\ t) = c a . (4) 

z^oo 

In (JH)- the field c is the cation concentration C, rescaled by a quantity R c that depends 
on the diffusivity and the mobility of the two species (cations and anions). Specifically, 
R c is the ratio between the ambipolar diffusion coefficient D and the constant t c , defined 
by 

\L c D a + p, a D c p c 
D = : tc = : , (5) 

He + Ha Hc + Ha 

where p a (p c ), D a (D c ) are the anion (cation) mobility and diffusivity, respectively. 
Other parameters in Eqs. (pQ) to (j3J are c° = R c C a exp[z c Fr]/(RT)}, c a = R c C a , 
T = 7/(i?T), Q = M/p = (1 — t c )/C a , and the mass transport coefficient hp = 
J exp[—bz c Fr]/(RT)]/(z c FC a R c ), where C a is the initial cation concentration, held 
fixed at infinity, M is the metal molar mass and p the aggregate mean density, T is 
temperature, 7 is the aggregate surface tension, R is the ideal gas constant and F 
is Faraday's constant, ez c is the cation charge, Jo is the exchange current density in 
equilibrium, and 77 is the overpotential from which a surface curvature contribution in 
the corresponding Butler- Volmer boundary condition [21 [251 EE] has been singled out, 
parameter b estimating the asymmetry of the energy barrier related to cation reduction. 

Eq. (JTJ) describes cation diffusion in the cell, with a noise q term that describes 
fluctuations in the (conserved) diffusive current. Eq. (T5]) is a mixed boundary condition 
at the aggregate surface ( that determines the efficiency at which the cation flux leads 
to actual aggregation. It originates from the standard Butler- Volmer condition, and 
allows for a noise term x accounting for fluctuations in attachment. In this equation 
d n represent the derivative along the normal component n of the interface. Eq. ([3]) 
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computes the growth velocity along the local normal direction V n = V • n, incorporating 
the contributions from the diffusive current mediated by the boundary condition (j2J), 
as well as from other diffusive currents at the surface, 3 S , in which fluctuations are 
described by the noise term p. The standard choice is that of surface diffusion [271 [28] . 
V s - J s = ^ 2 ^s-Ds7(^_b7')~ 1 V^C; where v s is the surface concentration of diffusing species 
and D s is their surface diffusivity. Here, V 2 denotes the surface Laplacian. Finally, 
the stochastic moving boundary problem is complete with the boundary condition (j3J) 
prescribing a fixed cation concentration at the anode. 

In the above model of ECD, in principle all parameters can be estimated in 
experiments. This is also the case for the noise terms: indeed, their amplitudes can 
be shown to be functions of the physical parameters above, once we take the noises 
as Gaussian distributed numbers that are uncorrelated in time and space, and make a 
local equilibrium approximation [10]. The mathematical problem posed by Eqs. (CQ) to 
PI is analytically very hard to work with. However, substantial progress can be made 
by means of a small slope, approximation |VC| «C 1 [291 ESS HD] • This approach allows 
to obtain a closed evolution equation for the Fourier modes (k(t) of disturbances of the 
surface height away from a flat interface. Generically, such an equation takes the form 



where is an additive noise term including both conserved and non-conserved 

contributions, and A/"[C]fc is the Fourier transform of jV[C] — {d x () 2 i namely, of 
the celebrated Kardar-Parisi-Zhang (KPZ) nonlinearity that is the landmark of non- 
conserved interface growth. Here, V is the mean velocity of a flat interface, while 
the dispersion relation u>k is a function of wave-vector magnitude k that provides the 
(linear) rate of growth or decay of a periodic perturbation Ck{t) of the flat solution. In 
the simplest long- wavelength limit, the dispersion relation becomes a polynomial in \k\, 
with two important limits, namely 



where B = Q 2, yi> s D s /(RT), and A = 1 — do/lo with do = TQ the capillary length 
and Id = D/V the diffusion length, the condition do <C Id applying in most physical 
cases. As we see from fl7J, in all cases the dispersion relation predicts a morphological 
instability in which there is a band of linearly unstable Fourier modes whose amplitudes 
grow exponentially in time. One of these modes grows at a faster rate than all the other, 
providing the selection of a characteristic length-scale in the surface and formation of a 
periodic pattern. However, some features of the instability change as a function of the 
mass transfer parameter ho- Thus, when attachment kinetics is very fast (ko oo) 
the dispersion relation has the classic shape of the Mullins-Sekerka (MS) instability. In 
this case, Eq. ([6j) is non-local in space, as a reflection of the strong implicit shadowing 
effects associated with diffusion. In the case of slow kinetics and, as long as do < Id, 



[221 una nig 



d t Ck(t) = WfcOtCt) + yAA[C] fe + Vk(t), 



(6) 




k D ^oo MS 
ko < oo KS 
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we still have a morphological instability, but of a local nature, Eq. ([2D becoming the 
celebrated (stochastic) Kuramoto-Sivashinsky (KS) equation in real space. 

Remarkably, the KPZ nonlinearity is able to tame the morphological instabilities 
in and induce, albeit after relatively long transients, kinetic roughening properties 
at long length scales pU]. These are those of the KPZ equation for slow attachment 
kinetics, while they correspond to a different universality class for kp — > oo. Thus, 
one would expect KPZ scaling for a wide range of parameters, which has indeed been 
found experimentally within ECD growth [31], but not with the generality that was 
expected [30], [32] . One possibility is that the crossovers induced by the instabilities are 
so long lived as to prevent actual asymptotic scaling from being accessed by experiments. 
Such a possibility seems to occur in conceptually similar growth techniques, such as 
Chemical Vapor Deposition [33] [34] . However, standard symmetry arguments [6] lead 
to the expectation of generic KPZ behavior for such non-conserved growth phenomena 
as provided by ECD systems. Thus, in order to understand this disagreement, in the 
next section we undertake a detailed comparison between the small slope approximation 
and ECD growth experiments. 



3. Dispersion Relation and Small Slope Approximation 

The simplest way to improve upon the analytical results of the previous section is to 
realize that Equation ([7]) merely provides a long-wavelength approximation of the exact 
linear dispersion relation for the moving boundary problem. We can consider the full 
shape of the dispersion relation, that turns out to be an implicit (rather than explicit) 
function of wave-vector, namely [35] 



1 + 



1 + iu k + k 2 



V2\ (_ V . d f2 



^ - = ^ + 4 ' ( 8 ) 



where k = 2Dk/V and Cj^. = ADu^/V 2 . Using this formula, we can see the extent to 
which surface kinetics modifies the functional form of the linear dispersion relation. 
Actually, several experimental studies of ECD have already tried to estimate the 
functional form of the dispersion relation in compact aggregates. An example taken 
from the literature is the growth of quasi-two-dimensional Cu and Ag branches by de 
Bruyn [36] [37]. In [37], the branches grow in the cathode due to the electrodeposition 
of ions created from the CUSO4 aqueous solution. The analysis of the early stages of 
the branch growth allows the author to fit the experimental dispersion relation to 

, l\k\(l-rk 2 ) 

Uk = 1 + sk ' (9) 
where q, r and s depend on the properties of the electrolyte and the deposited metal, 
the ion concentration, and the applied current. From the reported data and figures 
in Ref. [37J we can extract the information summarized in Table [H where we also 
include data from the experiments in [3B1 EH]- We need to supplement these values 
with reasonable estimates for the capillary length, <i and the mass transfer coefficient, 
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Parameter 


Ref. [37J 


Ref. [38J 


Ref. [39J 


C a [mol/cm 3 ] 


1 x 10" 4 


2 x 10~ 3 


5 x 10" 4 


I [mA] 


1.4225 


NR 


NR 


J [mA/cm 2 ] 


3 x 10~ 2 


3 x 10~ 2 


3 x 10 -2 


J [mA/cm 2 


11.16 


NR 


65.00 


D c [cm 2 /s 


0.720 x 10~ 5 


0.720 x 10" 5 


0.720 x 10~ 5 


D a [ 


cm 2 /s] 


1.065 x 10~ 5 


1.065 x 10" 5 


1.902 x 10~ 5 




2 


2 


2 


V 


cm/s 


3.45 x 10" 4 


« 1.00 x 10" 6 


4.87 x 10" 4 


Id 


cm 


2.49 x 10~ 2 


~ 8.59 


2.15 x 10~ 2 



Table 1. Experimental parameter values estimated from Refs. [37] (second column), 
[35] (third column), and [3J5] (4th. column). NR stands for not reported. 



kr>. We estimate do from the experiments by Kahanda et al. [38] (and summarized 
in Table [U see also below) and we rescale this value for the right concentration to 
do = 2 x 10~ 5 cm. Finally, kjj can be obtained from the location of the maximum in the 
dispersion relation. Thus we can estimate hp to be in the range (due to the experimental 
uncertainty) 4.44 — 5.56 x 10~ 4 cm/s. 




Figure 1. Linear dispersion relation, uj^, of Cu electrodeposition. Squares provide 
the experimental values taken from Ref. ,37] , and the red line is the fit to © made in 
this reference. The other lines are numerical solutions of ([5]) using the parameters in 
the second column in Table [1] for different values of kr> as indicated in the legend. 

These data are compatible with the assumptions made in the derivation of the 
moving boundary problem introduced in Ref. (TO]. For instance, the diffusion length 
is larger than the branch width (typically 0.1 mm in these experiments) but the mean 
growth velocity is comparable to the mass transport coefficient, and surface kinetics 
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become relevant. For this reason, the growth process is not diffusion limited and the 
dispersion relation is very different from the Mullins-Sekerka one. Similarly, the estimate 
for kjj implies that the applied overpotential rj varies between —0.160 V and —0.154 V, 
which are plausible enough. With this information we can now proceed to recast the 
experimental results in Ref. [37] as shown in Fig. [U The agreement between our theory 
(JEJ) and experiment is remarkable, taking into account that we are not performing any 
fit of the data. Moreover, the long wavelength approximations (J7J are not able to fit 
the experimental data. 

We have performed a numerical integration of an improved version of Eq. ([6]) in 
which the full linear dispersion relation (jBJ) is employed, using a stochastic pseudospectral 
scheme, see e.g. [^Ql HQ] - Results are shown in Fig.[2]for the time evolution of the surface 
roughness W(t) (root mean square deviation) and the surface structure factor or power 
spectral density (PSD), S(k,t) = (|Cfc(^)| 2 ), of the surface height, ((x,t). As we see, the 
roughness increases very fast (close to W(t) ~ t 18 ) within an intermediate time regime, 
after which saturation to a stationary value is achieved. 




t[s] 



Figure 2. Time evolution of the surface roughness and PSD (insets) using $Q with 
u>k given by ([5]) (left panel) and (TUT)) for [i — 0.75 (right panel). 

As for the time evolution of the PSD curves, they follow a standard Family- 
Vicsek (FV) behaviour [6], leading to a global roughness exponent a = 1.25, from the 
asymptotic behaviour S(k) ~ k~( 2a+d ^ for a ID (d = 1) interface. The time evolution of 
the interface height in real space can be appreciated in movie 1 in the supplementary 
material. As is clear from the images, an initially disordered interface develops mound- 
like structures that coarsen with time, leading to a long time super-rough morphology. 
This is inconsistent with the above expectation of KPZ scaling for the present case since, 
taking into account the values of V and ko, a KS behavior is expected for (IE]). In the 
inset of Fig. [3] we show the full dispersion relation ([8]) for sufficiently small values of 
wave-vector at which behaviour is indeed of the KS type. This apparent contradiction 
can be explained taking into account that such a KS shape (that holds for length scales 
in the cm range) cannot be observed in practice for the physical length-scales that 
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Figure 3. Full linear dispersion relation as shown in Fig.[T]for hp = 5 • 10~ 4 cm s _1 . 
For the meaning of the solid line see the main text. The inset is a zoom showing the 
Wfe ~ fc 2 behaviour for sufficiently small k values. 



are accessible in the experiment. As seen in the figure, the effective behavior of the 
dispersion relation is not of the KS type but, rather, ujk behaves (sub)linearly for the 
smallest accesible k values. This is made explicit in the main panel of Fig. [3j where a 
straight line is shown for reference at the smallest wave-vectors. Thus, for all practical 
purposes the linear dispersion relation that is experimentally operative does not behave 
as ~ k 2 but, rather, is closer to a form such as 

u ll (k)=b li \k\' t -b 2 k 2 , (10) 

where and 62 are positive, effective parameters, and < /i < 1. 

Taking into account the previous paragraphs, we see that ( flOl) summarizes the 
main features of the linear dispersion relation for many experimental systems, namely, 
(i) fast growth at the smallest available k values, ujk ~ {k^, (ii) single finite maximum 
at intermediate k values, and (iii) sufficiently fast decay at large wave-vectors, at least 
of the form Uk ~ —k 2 or faster (i.e. with a larger exponent value). Thus, a simple way to 
explore the dynamics as predicted by the small slope approximation with an improved 
dispersion relation consists in considering Eq. ([6]), but using ffTOl) rather than ([7j). We 
have carried out such type of study for several values of [l. For the sake of reference, 
movie 2 in the supplementary material shows the evolution of an interface obeying the 
KS equation, while movie 3 displays the analogous evolution but as described by Eq. 
(EJ) using the effective dispersion relation ( flOl) for \i = 0.75. It is apparent that the 
latter resembles much more closely the dynamics associated with the full dispersion 
relation (jHJ). Quantitative support is provided by the right panel of Fig. El where the 
evolution of the roughness and PSD is shown for Eq. with (TTDT) for \i = 0.75. In 
the inset, we can appreciate the formation of a periodic surface (signalled as a peak 
in the PSD) at initial times that grows very rapidly in amplitude, as confirmed by the 
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exponential evolution of the global surface roughness W(t), shown in the main panel 
of Fig. [21 for times t < 15. For long enough times, this peak smears out, once the 
lateral correlation length is larger than the pattern wavelength. That marks the onset 
for nonlinear effects that induce power-law behavior (kinetic roughening) both in the 
roughness and in the PSD, approximately as W(t) ~ i 1,85 and S(k) ~ l//c 3,5 , wherefrom 
we conclude that the associated critical exponents have values f3 = 1.85 and a = 1.25 
(thus, z = a /(3 = 0.68, not far from z = \x = 0.75, an exponent relation which at 
least holds for ( Tl0|) when fi = 1 [TU].) We find it remarkable that the KPZ nonlinearity 
is able to stabilize the initial morphological instability of Eq. ([6]) even in the case of 
a dispersion relation such as 0^=0.75 that induces such a fast rate of growth. Notice 
z = 0.68 < 1 corresponds effectively to correlations building up at a rate that is even 
faster than ballistic. Moreover, the height profile develops values of the slope at long 
times that are no longer compatible with the |VC(£)| <C 1 approximation that led to 
the derivation of Eq. ([6]) in the first place. This requires us to reconsider the original 
moving boundary problem from a different perspective. 

4. Multivalued interfaces 

In order to go beyond the small-slope approximation one has to face the numerical 
integration of the full moving boundary problem (pQ)- <^j. As mentioned in the 
introduction, a method of choice to this end is the so-called phase-field or diffuse interface 
approach [11]. In this method, one introduces an auxiliary (phase) field <fi(r,t) that is 
coupled to the physical concentration field c(r, t) in such a way as to reproduce the 
original moving boundary problem asymptotically, within a perturbation expansion in 
the small parameter e = W/Iq. Here, W denotes the typical size of the region within 
which the value of <fi(r,t) changes between <fi = — 1 in the vapor phase and <fi = +1 in 
the solid phase. The level set = will thus provide the locus of the aggregate surface, 
avoiding the need of explicit front tracking. For our present physical problem, such an 
asymptotic equivalence has been recently established [35] for the following phase-field 
reformulation of the stochastic system (CQ)-(J4]), 

rd t( f> = W 2 V 2 (P - d4f(4>) - Xug(4>)} + (11) 

d t u = V- {Dq{4>) Wu - j at ) - ~d t /i(0) - V • <t u , (12) 

where, u = fi(c — c° eq ), and cr u are noises whose variances are related to those of q 
and x i n ©-© [35] , and f(4>), g(4>), q{<j>), and h(<p) are auxiliary functions that are 
conveniently chosen (in order to achieve asymptotic convergence) as /(</>) = 4 /4 — 2 /2, 
g(4>) = <p - 20 3 /3 + 5 /5, q{4>) = (1 - 0)/2, and h{4>) = <p [41J. Moreover, the current 
j a t generates a so-called anti-trapping flux that is needed for technical reasons related 
to convergence [HI [35], and can be taken as j a t = 8^ 1 ^ 2 W(d t (p) V0/|V0|. Given model 
(JTT|)-(fr2l). one can finally prove [35] its convergence to the original moving boundary 
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0.5 mm 






Figure 4. Comparison of morphologies from the experiments on CUSO4 in Ref. [38] 
(left column, taken from the original paper by permission; reference bars represent 
0.5 mm) and phase-field simulations using parameters in Table [T] (right column). 
Top row corresponds to fast kinetics (77 = —0.50 V), and bottom row corresponds 
to intermediate kinetics {q = —0.35 V). 



problem in the thin interface limit e< 1, provided model parameters are related as 

A - W 
"0 — 

( T W\ 

k» =^[xW~ a2 D)' (13) 
where a± and a 2 are constants associated with the phase-field model, with values 
ai = 0.8839, 

a 2 = 0.6267. (14) 

There have been previous proposals of phase-field formulations of diffusive [121 1121 
HI] (and ballistic [HI HI]) growth systems of the type of our moving boundary problem. 
However, in these works the phase-field equations are proposed on a phenomenological 
basis, rather than being connected with a physical moving boundary problem as in our 
case. Thus, quantitative comparison is harder to achieve, although many qualitative 
morphological features can indeed be retrieved. 

We have performed numerical simulations of model (flTT) - ffT2l using physical 
parameter values that correspond to the experiments of Cu ECD by Kahanda et al. [38] 
(see parameters in Table [T|), for different surface kinetics conditions. Representative 
morphologies are shown in Fig. H] where the overpotential value is employed to tune the 
value of the kinetic coefficient hp- Notice that in Ref. [38] the average growth velocity 
is overestimated, having been measured in the nonlinear regime, which also overshoots 
the estimated value of the capillary length. For this reason, the value of the capillary 
length employed in our simulations (and provided in Table [1]) is estimated as smaller 
than that in [38]. As we can see, for the faster kinetics case (upper row in Fig. Hj), 
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indeed the interface becomes multivalued already at early times, leading to a relatively 
open, branched morphology As the overpotential increases (thereby decreasing the 
value of hp, lower row), the morphology is still multivalued, but more compact. We 
would like to stress the large degree of quantitative agreement between the simulated 
and experimental morphologies. 

Our phase-field model (and the corresponding moving boundary problem that it 
represents) has moreover a large degree of universality, in the sense that it can also 
reproduce quantitatively other ECD systems. An example is provided by Fig. El which 
shows simulations for ECD growth from a Cu(N0 3 ) 2 solution, as in Ref. [39], see 
experimental parameters in Table [TJ This dense branching morphology corresponds 
to a fast kinetics condition for which the single valued approximation breaks down 
already for short times. Fig. [6] collects the time evolution of the surface roughness and 
of the PSD function for the same parameters as in Fig. El The roughness grows rapidly 
within a transient regime, that is later followed by stabilization and rapid saturation 
to a stationary value. This result is similar to those obtained within the small slope 
approximation. However, the behavior of the PSD differs dramatically. Thus, for very 
small times the linear instability shows up in the finite maximum of the curves. For 
intermediate times the PSD's shift upwards with time, a behaviour that is associated 
with anomalous kinetic roughening properties [8j 06] that contrasts with the standard 
FV scaling found above in which the PSD values are time independent for length scales 
that are below the correlation length. Note that at these intermediate times the single 
valued approximation does not describe the features of the interface at all. Actually, 
the single valued approximation to the height field (not shown) presents abundant large 
jumps that are possibly responsible for the effective anomalous scaling [47] . However, 
for sufficiently long times branches grow laterally closing the interbranch spacing, and 
the validity of the single valued approximation is restored. Thus, the PSD curves for the 
longest times again display Family- Vicsek behavior for the largest k values. A crossover 
between a = 0.5 for small scales and a = (log) for intermediate scales can be read off 
from the figure. 

Note that, for the longest distances in our system, scale invariance seems to break 
down, since no clear power law can be distinguished for, say, k < 30 mm -1 ~ l D x . In 
order to properly characterize the behavior of the interface at these scales, simulations 
are required for system sizes that are much larger the diffusion length. We hope to 
report on this in the future. 

5. Discussion and Conclusions 

In this paper we have studied the physics of non-conserved growth from a diffusive 
source. We have focused on, both, morphological features and dynamical properties, 
making special emphasis in the emergence of a long time structure arising from a well 
characterized linearly unstable regime. Thus, we have dealt with the problem from two 
different but complementary points of view: the linearized dispersion relation obtained 
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Figure 5. Simulated morphology using parameters as in the experiment of Ref. |39j 
on Cu(NC>3)2, corresponding to a fast kinetics condition. Lower panel: zoom of a 
portion of the upper panel. 

from the moving boundary description of the problem (as presented in Ref. [10]) and 
the use of a phase-field equivalent description of the full original physical system. 

This complementary study has allowed us to understand the evolution of the 
problems in two well separated regimes: one at short times, well characterized by the 
linear dispersion relation, and another one at long times in which steep slopes may 
eventually develop and break down the simplistic description through a single valued 
height field. We want to emphasize that the comparison has been done quantitatively. 
Thus, in Sec. II we have shown how our theoretical predictions allow to explain some 
experimental results reported in [37] . At variance with the data analysis performed in 
this reference, our comparison is not based on a multidimensional fit of the parameters 
but arises, rather, from a prediction for the form of the dispersion relation compatible 



Kinetic 



14 




10" 



10" 10" 



io u 
t[s] 



10 10 



10 



l(f 

k [1/mm] 



10' 



10 



io 4 



10 



10 



10" 



10" 



o »3 



Figure 6. Numerical time evolution of the surface roughness (left panel) and PSD 
(right panel) for the same parameters as in Fig. [5l The PSD curves are for increasing 
times from bottom to top. Notice the time shift of the PSD curves for high values 
of k (small length scales) signalling anomalous scaling, see the text for details. For 
reference, the dashed line has slope —2 while the dash-dotted line has slope —1. 



with the available experimental data. 

In order to develop insight into the dynamics decribed by the complex full dispersion 
relation (J5J), we have proposed a simple description of the evolving interface in terms 
of an effective dispersion relation, in which the exponent \x of the destabilizing term 
comprises all the relevant information on the behavior of the morphological instability 
at the largest length scales accesible in the system. Hence, the simple description of 
the data given by Eq. (jHJ) using (jSJ) provides a promising testbed for future theoretical 
analysis of many problems sharing similar ingredients with ours. Another benefit of 
this characterization of the linear regime is to provide an explanation for the scarcity 
of experimental evidence of the celebrated KPZ scaling (or even its generalization for 
unstable systems given by the KS equation). As we have demonstrated, the difficulty 
of actually observing the experimental systems in the KPZ universality class is due 
to the fact that the physical conditions for the observation of that behavior cannot 
be fulfilled for realistic (experimentally motivated) values of the parameters, as they 
typically require system sizes that are beyond experimental feasibility. 

Besides, with the help of the phase-field model we have also performed a systematic 
analysis of the experiments reported in [38]. Our comparison not only reproduces 
qualitatively well the morphologies obtained in those experiments but also predicts 
the existence of a scaling regime characterized by anomalous scaling. This anomalous 
scaling deserves further analysis but, preliminarily, can be understood in terms of the 
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emergence of large slopes for times following appearance of overhangs and multivalued 
interfaces. As we have shown, our model sheds some light on the characterization of the 
experiments reported in Refs. [381 139] and the role of the different mechanisms involved. 
The quantitative estimation of the parameters appearing in our theoretical framework 
leads also to remarkable agreement with the morphological structure of the observed 
deposits. 

Finally, we want to stress that, although the formulation of the electrochemical 
deposition process with a phase-field theory has been addressed before in the literature 
[121 H3j HU H5], our model is not constructed phenomenologically, but is, rather, an 
equivalent formulation of the physical equations (diffusion, electrochemical reduction 
at the cathode, etc.) that reduces to them in the thin-interface limit. In addition, 
our formulation allows us to do a quantitative comparison with experiments which is 
far from the capabilities of the mentioned theories. Notwithstanding, the promising 
results presented above still need more detailed comparison and analysis (through 
more extensive numerical experiments), which will be the aim of future work. This 
will focus particularly on the appearance of the anomalous scaling regime and, from 
the morphological point of view, on the influence, cooperation and/or competition of 
different physical mechanisms that are naturally included in our theory. 
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